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فناوری‌های نوین در مهندسی برق و سیستم انرژی سبز 
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یک طرح بهینه جهت هماهنگی بازیابی بار و کاهش زاوبه فاز ایستا با استفاده از 
الگوریتم ۲1,80" 


هادی حسین پور ‏ دانشجوی دکتری» محمدرضا اسماعیلی " دکتری ۰ امین خدابخشیان ‏ استاد 


۱- دانشکده مهندسی برق, دانشگاه اصفهان, اصفهان ایران 
۲- شرکت برق منطفه‌ای اصفهان, اصفهان, ایران 


چکیده: علیرغم وجوه کلیه تجهیزات کنترلی و حفاتی ذر یک سیستم قمرت, احتمال وقوع یک خاموشی امری اجثتاب ناپذیر انست. 
بنابراین. فرآیند بازیابی یکی از مهمترین دغدغه‌های بهره برداران سیستم است تا بتوانند در کمترین زمان ممکن خسارات ناشی از آن را 
کاهش دهند. در فرآیند بازیابی موازی ابتدا جزایر مورد نظر تشکیل شده و سپس بار هر جزیره به طور جداگانه و بصورت همزمان بازیابی 
می‌شود. در مرحله بعد. جزایر تشکیل شده با رعایت حداقل مقدار زاوبه فاز ایستا با یکدیگر سنکرون شوند. برای انجام اين کار. یک طرح 
چند هدفه بهینه در اين مقاله تعریف شده است تا مسائل بازیابی بار و کاهش ٩۳۸‏ را به طور هماهنگ بهینه سازی نماید. توابع هدف 
مدل پیشنهادی شامل به حداقل رساندن زاویه فاز ایستا و به حداقل رساندن انرژی تامین نشده است که با در نظر گرفتن محدودیت‌های 
مورد نظر بهینه می‌شوند. در اين راستا از الگوریتم بهینه سازی آموزش و یادگیری (11,20) به عنوان تکنیک پیشنهادی استفاده شده 
و با برخی از الگوریتم‌های هوشمند مقایسه شده است. شبیه‌سازی‌ها با ایجاد ارتباط بین دو نرم‌افزار ]۷۸1 و 1109511/۳۲ 
انجام می‌شود. نتایج به‌دست آمده نشان‌دهنده کارایی مدل پیشنهادی برای دستیابی به اهداف ذکر شده است. 


واژه های کلیدی: بازیابی سیستم قدرت. زاویه فاز ایستاء بازیابی بار. انرژی تامین نشده 


تاریخ ارسال مقاله: ۱۴۰۰/۱۲/۲۹ 

تاریخ بازنگری مقاله: ۱۴۰۱/۰۳/۰۸ 

تاریخ پذیرش مقاله: ۱۴۰۱/۰۳/۲۵ 

نویسنده‌ی مسئول: دکتر محمدرضا اسماعیلی 13606210211.000 1513012111 
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هادی حسین‌پوره محمدرضا اسماعیلی امین خدابخشیان 


۱- مقدمه 


بازیابی سیستم قدرت پس از بروز خاموشی سراسری يا جزئی یکی از مسائل مهم شبکه‌های قدرت می‌باشد. شناسایی و انجام پیکربندی 
مجدد شبکه به‌صورت اصولی باعث می‌شود تا شبکه به سرعت و به بهترین صورت ممکن بازیابی شده و مدت زمان خاموشی برای 
مشترکین به حداقل برسد [۱]. از آن‌جایی که وابستگی جوامع مدرن به انرژی الکتریکی به سرعت در حال افزايش می‌باشد تقاضا برای 
دریافت مطمئن این انرژی حیاتی در حال افزايش است که این امر باعث پیچیده‌تر شدن عملکرد سیستم‌های قدرت و احتمال بروز خطا 
و خاموشی در شبکه می‌شود. بروز یک خاموشی سراسری در شبکه برق می‌تواند خسارات جبران ناپذیری به یک جامعه تحمیل نماید. از 
این‌رو یکی از مهمترین دغدغه های بهره برداران سیستم قدرت برگرداندن وضعیت سیستم به حالت عادی یا همان بازیابی سیستم قدرت 
می‌باشد [۲]. هدف نهایی در فرآیند بازیابی سیستم قدرت. بازیابی بار در کمترین زمان ممکن و کاهش خسارت‌های مربوط به مصرف 
کننده سیستم قدرت می‌باشد [۲]. تحقیقات متعددی در خصوص طراحی برنامه‌های بازیابی سیستم قدرت در سه بخش خودراه اندازی» 
پیکربندی و بازیابی بار توسط محققین ارائه گردیده است [۴-۶]. در مرجع [۴] یک روش. مبتنی بر تئوری گراف ارائه شده است که در 
آن سیستم قدرت به‌صورت یک گراف در نظر گرفته شده و مسیرهای بازیابی از روی این گراف انتخاب می‌گردد. در مرجع [۵] یک مدل 
بهینه‌سازی ترکیبی با در نظر گرفتن مسئله تعیین توالی حذف بار در آخرین مرحله از پیکربندی مجدد شبکه به کار گرفته شده است. 
در مرجع [۶] یک روش برای تخمین زمان بازیابی در تمامی فازهای بازیابی بر اساس تکنیک یاه گیری ماشین ارائه شده است: 

مرحله بازیابی بار در بخش‌های انتهایی برنامه بازیابی حساس‌ترین مرحله به شمار می رود. در این مرحله بهره برداران با دو چالش 
بزرگ شامل مساله کاهش انرژی تامین نشده و کاهش زاویه فاز ایستا در زمان سنکرون سازی جزایر روبرو هستند. مساله کاهش زاویه 
فاز ایستا در زمان سنکرون سازی دو زیرسیستم یکی از پیچیده‌ترین اقدامات در بهره برداری سیستم قدرت می‌باشد که میزان تولید و 
بار سیستم در نقاط مختلف شبکه به طور قابل ملاحظه‌ای بر روی آن تاثیر گذار است [۷-۱۹]. لذا همواره تدابیری اتخاذ می‌شود تا زاویه 
فاز ایستا نیز کاهش یابد و عمل اتصال خطوط میان نواحی و یا باس ها امکان‌پذیر باشد [۷ و ۸]. اختلاف راویه فاز ولتاژ بین دو باس در 
سیستم قدرت به عوامل مختلفی از جمله. ولتاژ نامی. حالت عملکرد شبکه و محل قرارگیری کلید قدرت مورد نظر بستگی دارد .]٩[‏ 
معمولا این زاویه تحت شرایط بازیابی سیستم قدرت باید به مقدار ۲۰ درجه یا کم‌تر کاهش یابد. یکی از راهکارهای کاربردی کاهش 
اختلاف زاویه فاز ایستاه تغییر در مقدار توان تولیدی ژنراتورهای سیستم قدرت است [۱۰]. مسئله کاهش ٩۳۸‏ در مرجع [۱۱] از طریق 


8 یک طرح بهینه جهت هماهنگی بازیابی بار و کاهش زاویه فاز ایستا با استفاده از الگوریتم (11,۳6[" 


برنامه ریزی تولید نیروگاه‌ها بدون نیاز به تغییر در بار حل شده است. در [۱۲ به‌منظور کاهش زاوبه فاز ایستا از یک استراتژی مبتنی 
بر خطی‌سازی و برنامه‌ریزی درجه دو استفاده شده که در آن توان تولیدی واحدها. پخش بار و محدودیت‌های زاوبه فاز ایستا در نظر 
گرفته شده است. در[۱۳] روش مبتنی بر برنامه‌ریزی مجدد توان تولیدی سیستم قدرت به‌منظور کاهش زاویه فاز ایستا استفاده شده 
است. در [۱۴]» برای دستیابی به زاویه فاز ایستا مطلوب. ولتاز ترمینال ژنراتور به‌عنوان متغیرهای کنترل در نظر گرفته شده است. در 
[۱۵] به بررسی تاثیر ادوات 3۸015 در کاهش زاویه فاز ایستا پرداخته شده است. مرجع [۱۶] امکان استفاده از ۱۳۳ را برای کاهش 
زاوبه فاز ایستا در سیستم قدرت مورد بررسی قرار داده است. یکی از تجهیزات قابل کنترل در سیستم‌های قدرت الکتریکی که برای 
تغییر اختلاف زاویه فاز ایستا مورد استفاده قرار می‌گیرد» تپ ترانسفورماتورها است. با تغییر در ولتاژ شین‌های شبکه منجر به ایجاد 
تغییرات در زاوبه فاز شین‌ها می‌گردد [۱۷]. در بسیاری از شرایط تغییر بیش از حد توان تولیدی زنراتورهای شبکه در جهت دستیابی به 
حد مجاز اختلاف زاویه فاز ایستا فراهم نیست. به همین منظور پيشنهاد استفاده از برنامه حذف بار به‌عنوان راهکاری برای حل مشکل 
بیان شده مطرح می‌شود [۱۷]. در مرجع [۱۸] یک روش جدید برای کاهش ۹۳۸" ارائه شده است که در آن برداشت بار و افزايش توان 
تولیدی به‌عنوان یک ابزار کنترلی در سیستم در نظر گرفته شده است. 

در تحقیقات صورت گرفته مساله کاهش زاویه فاز ایستا بدون در نظر گرفتن مساله کاهش انرژی تامین نشده در حین فرایند بازیابی 
بار مورد بررسی قرار گرفته است. در این مقاله با ارائه یک طرح چند هدفه دو مساله کاهش زاویه فاز ایستا و کاهش آنرژی تامین نشده 
به طور همزمان در حین فرآیند بازیابی بار بهینه‌سازی می گردد. نوآوری های این تحقیق عبارتند از: 

۶ هماهنگی دو مساله کاهش زاوبه فاز ایستا و کاهش انرژی تامین نشده در حین فرآیند بازیابی بار. 

9 استفاده از روش هوشمند 1120" جهت بهینه سازی مساله پیشنهادی. 
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۲- فرمول‌بندی مساله پیشنهادی 

از مجموعه فرآیند بازیابی سیستم قدرت. بازیابی بار با هدف کاهش زاویه فاز ایستا و انرژی تامین نشده بصورت یک مسئله بهینه‌سازی 
مورد بررسی قرار می‌گیرد. میزان انرژی تامین نشده و ماکزیمم زاویه فاز ایستا به‌عنوان اهداف بهینه‌سازی در هر دو مسئله در نظر گرفته 
انا مه ودیت‌های خربوط به کوان‌های اخشویی را وبا ها وان قالیتی تا د ها پیره نظ نار او آسای کهقر باربای 
شبکه متغیرها زیاد می‌باشند. جهت بهینه‌سازی از روش‌های فرا ابتکاری و روش‌های هوشمند استفاده می‌کنند [۲۰]. در اين پژوهش. 
الگوریتم بهینه‌سازی مبتنی بر آموزش و یادگیری (11.1300) جهت حل مسئله بهینه‌سازی پيشنهاد شده است و از الگوریتم‌های اجتماع 
ذرات((۳۹6). ژنتیک(7۸)) و الگوریتم کرم شب‌تاب (۳۸) جهت مقایسه با الگوریتم پیشنهادی استفاده شده است. با توجه به پیچیددگی 
شبکه قدرت. بهینه‌سازی باید در نرم افزار متلب انجام شود و همچنین جهت جزیره سازی شبکه به‌منظور بازیابی بار و انجام پخش بار 
تالم انقر دق اقا میکتالانت انشاهه کقد تا امتان حاصل نموه اسش سک آموه تیک ری اسف به ال ای شک 
دارد. براین اساس در این مقاله. با برقراری ارتباط بین این دو نرم افزار عمل بازیابی صورت گرفته است. در ادامه دو تابع هدف مجموع 
انرژی تامین نشده و ماکزیمم زاویه فاز ایستا معرفی و فرمول بندی آن‌ها ارائه می گردد. 

کاهش زاوبه فاز ایستا در زمان سنکرون‌سازی دو زیرسیستم به‌عنوان یکی از مهم‌ترین معیارها تلقی می‌شود. یکی از مهم‌ترین عوامل 
تاثیرگذار بر روی این معیار مقدار توان اکتیو تزریق شده با مصرف شده در هر باس‌بار می‌باشد. برای بررسی این موضوع. ساختار 
ساده‌ای از یک خط انتقال در شکل (۱) نشان داده شده است که در آن 2 به‌عنوان امپدانس خط انتقال در نظر گرفته شده و از 
ادمتاشسن بط انا صقن شده استقر ایخ شک ولفاه انتهای خط اقا پمتوای فادیی مرجم(0ا او ملاه ایتدای خظ 
با اختلاف یک زاویه فاز 0 (1/:20-و/)) نسبت به ولتاژ انتهای خط بصورت یک فازور پیش‌فاز در نظر گرفته شده است. توان ظاهری 
ارسالی از ابتدای خط و توان ظاهری دریافتی در انتهای خط به‌ترتیب بصورت ,60 + ,۲ < ,۵ و ,40 + ,]2 , مفروض می‌باشند. 
بر این اساس. جریان عبوری از خط انتقال باعث افت ولتاژ در مسیر می‌شود که رابطه آن بر اساس مقادیر توان و ولتاژ سمت انتهای 
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خط به‌صورت معادله (۱) نوشته می‌شود [۲۱]. 


۳1-۲ 0+ 
)۱( بلاط + ب را ۳ در کت 


۳ ۳ 


7:2 


فاز ولتاز انتهای خط انتقال به‌صورت رابطه (۲) نوشته می‌شود: 


۸۹0 

۹3 حسفتت-- ما۸ < 0 , لا +( ,۵+ 00 ,0 
نا خاط ‏ با ۱ 

مقدار مقاومت در معادل (۱) صرف‌نظر نمود و مقادیر افت ولتاژ حقیقی و موهومی را به‌صورت زیر بازنویسی کرد: 


۳۹ 0 
۸0 ۳ 
0 0] 


1 
/ ۲ 


با توجه به روابط (۲) و (۲) مقدار زاویه فاز ابتدای خط انتقال (در سمت ارسال کننده توان) مطابق با معادله (۴) به‌دست می‌آید: 


۲ 
)۳( سس ۸1۵0 - 0 


۱ د 1 

و 
با فرض انجام کنترل ولتاژ در شبکه و ثابت ماندن توان راکتیو تزریقی در هر باس تغییرات زاوبه فاز ( ۸0) به‌طور مستقیم متناسب با 
تغییرات توان اکتیو خواهد بود. مقدار تغییرات زاویه فاز در هر باس‌بار با توجه به معادله (۴) به شکل زیر خواهد بود: 


دهمکمی 

مکح 

مجججح 
سس ۰ ۲ ۰ اه " ۲۷ اک 
فناوری‌های نوین در مهندسی برق و سیستم آنرژی سبز. سال اول» شماره ۰۱ بهار ۱۴۰۱ شحییع 


(۵) مرس --۸0 


شکل(۱): مدار معادل یک خط انتقال 


با توجه به رابطه (۵) می‌توان بیان نمود که یکی از موثرترین ابزارهای تغییر زاویه فاز ایستا مابین دو زیرسیستم» میزان بار بازیابی شده 
بر روی سیستم در هر مرحله از فرایند باررگذاری می‌باشد. بر این اساس اولین تابع هدف بصورت حداقل‌سازی بیشترین زاویه فاز ایستا 
(5۳۸) مابین زیرسیستم‌های تشکیل شده و قبل از عملیات سنکرون‌سازی زیرسیستم‌ها تعریف می‌شود. در اين راستاء ٩1۸‏ بین دو 
باس انتهایی زیرسیستم‌های مجاور پس از اجرای عملیات پخش بار در هر دوره از فرآیند بازیابی بار به‌دست می‌آید. برای حل این 
مساله» ابتدا بیشترین مقدار ٩۳۸‏ در بین تمام خطوط ارتباطی بین زیرسیستم‌ها محاسبه میگردد. در نهایت» حداقل‌سازی این تابع 
هدف با کمک یک الگوریتم اجرا شده تا متغیرهای تصمیم گیری (مقدار بهینه بار قابل بازیابی به‌صورت بهینه به‌دست آیند. فرمول‌بندی 
تابع هدف به‌صورت رابطه (۶) نوشته می‌شود: 


8 یک طرح بهینه جهت هماهنگی بازیابی بار و کاهش زاویه فاز ایستا با استفاده از الگوریتم (11,۳6[" 


(۶ : 7 ۱9 
اه ور و ا-- 


فص رگ 
وه ود 4 ی رگ صنصه 14 

-1۳9* 
6, > ۵ 


صنصدط برگ 


حقص 4 


که در معادله (۶» پارامترهای ۰ زگ ریک رگ رگ ره 


۱ ۳ ۳ به‌ترتیب بیانگر خط انتقال ام بین دو 


5۳ من رگ 14 ها منم وا ٩‏ 
باس انتهایی و /. زاویه فاز ایستا بین دو باس انتهایی 7 و # توان تولیدی ژنراتور 7 ام. حداکثر توان تولیدی زنراتور 77ام. حداقل توان 
زاویه فاز ایستا بین دو باس انتهایی 7و می‌باشند. 


۱-۲- حداقل‌سازی مجموع انرژی تامین نشده (1۵۷5) 
دومین تابع هدف. حداقل‌سازی بیشترین مقدار مربوط به ٩‏ می‌باشد. که با استفاده از الگوریتم ابتکاری بهینه‌سازی می‌شود تا 
متغیرهای تصمیم گیری (مقدار بهینه بار قابل بازیابی) به‌صورت بهینه به‌دست آیند. معادله مربوط به این تابع هدف در رابطه زیر آورده 


شده است: 
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[(۸/۵)۱5 ی 7 صنصر 
ط 
)۳( زگ گ سم رگ وییز1 
او یا ر یی نانوی 
۷وو2ر1ع< ,ی ۲ > 


قص رگ منک 
رماع > وک ط 
هه 0 «ن ۵ 


ی ۳۹ 5 


در رابطه (0۷, مقادیر مربوط به 5( به‌صورت یک معادله در رابطه (۸) آورده شده است [۲۱]. 


۲۷5 < 


بص| 


1 
)+( 6 رت ارم 


1 <1 


1 


در رابطه (۸) پارامتر آ ۰0 "رگ و #به‌ترتيب نشان‌دهنده کل زمان بازیابی» کل بار موجود در سیستم. ماکزیمم بار قابل بازیابی و 
بار بازیابی شده در زمان نا می‌باشند. 


۲-۲- الگوربتم پیشنهادی 

الگوریتم بهینه‌سازی برپایه آموزش و یادگیری. مانند الگوریتم‌های بهینه‌سازی ۸ ۲30 یک روش مبتنی بر جمعیت است که 
گروهی از راه حل ها را برای ادامه فرآیند بهینه پیاده سازی می‌کند. بسیاری از روش‌های بهینه‌سازی نیاز به پارامترهای الگوریتم دارند 
کی که نا ترچ کار گنه بارس های ایا ای ترش ین بر انار ار ۳۱0 یه رای بافگیری؛: 
تغییر وزن و حداکثر مقدار سرعت نیاز دارد. برخلاف سایر تکنیک‌های بهینه‌سازی. 11,20 به تنظیم پارامترهای الگوریتمی نیاز ندارد. 
بنابراین اجرای ساده‌تر می‌شود. ارآ[ مانند ۳۵۵ از بهترین راه حل تکرار برای تغییر راه حل موجود در جمعیت استفاده می کند. 
در نتیجه نرخ همگرایی را افزايش می‌دهد. 1126 جمعیت را مانند ۸261 تقسیم نمی کند. همانند 0۸ که از انتخاب. ادغام و جهش 
استفاده 9 0 که از زنبورهای شاغل. تماشاگر 9 پیش‌اهنگی استفاده می کند. ]1 از دو مرحله مختلف معلم 9 دانش آموز 9 
همچنین از مقدار متوسط جمعبت برای به روزرسانی راه حل اسستفاده می کند. در نتیحه نرخ همگرایی ر افزایش می‌دهد. 0 ]1 
جمعبت ر مانند 0 تقسیم نمی کند. همانند ۵ که از انتخاب. ادغام و9 جهش استفاده 9 0 که از زنبورهای شاغل. تماشاگر و9 
پیش‌اهنگی استفاده می کند. 0ص از دو مرحله مختلف معلم 9 دانش آموز و همچنین از مقدار متوسط جمعبت برای به روزرسانی 
راه حل استفاده می کند. در این پژوهش از الگوریتم دص آبه‌عنوان الگوریتم پیشنهادی برای بهینه‌سازی استفاده شده است. شکل (۲( 
داده ها؛ ۳( تعیین حدود متغیرها؛ ۳( فرمول‌بندی مسئله 9 محاسبه توابع هدف‌های بیش‌ترین زاوبه فاز ایستا مابین زیرسیستم‌ها 9 
فختیون آترژی امین تفه ۴) نحل مشاه پرتامه ریتعنبانیای با ۵) کت کرفی سصنوذیت‌های سبعله ۶ تتایه خر وج 


۳-۲- برقراری ارتباط بین نرم افزار متلب و دیگسایلنت 


حل مسائل بهینه‌سازی در شبکه‌های قدرت. به دلیل گستردگی شبکه‌ها و پیاده‌سازی آن‌ها در محیط‌های پیچیده نرم افزاری» همواره 
یکی از دغدغه‌ها و چالش‌های پیش‌روی بهره‌برداران سیستم‌های قدرت بوده است. اگر چه انتقال اطلاعات از محیط نرم افزارهایی مانند 
۳ ]1281 به محیط نرم افزاری مانند 11,۸۳ 1۷۸ نیازمند تلاش و دقت فراوان است. اما از آن جا که بهینه سازی باید در محیط 
1 انجام شود. بهترین راه مورد استفاده توسط محققان در زمان انجام مطالعات شبکه . انجام پخش بار در محیط نرم افزار 
دیگسایلنت و بهینه‌سازی توسط نرم افزار متلب می‌باشد. در این پژوهش سعی شده است که از امکانات نرم افزارهای ]112911۳۱ و 
به‌طی. کایل اسفایه کزهوغا پایه‌های پشدسه آمهه شین امش یه حال اف که پاش در این ریش عل 


ارتباط بین هو ترم‌افزاز اسفاده شده است, مطالیات. بخ بار خر مخیط نرم افزا ۲۲ ااقعاط و انجام مراحل پپیته‌ساری هر محیط 
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نرم افزار 1۷۸11۸7 انجام می‌گیرد. در این شرایط نیازمند هیچ گونه ساده سازی در مورد شبکه مورد مطالعه نیست و تمامی حالت‌هاء 
ناشی از حالت واقعی شبکه هستند. محیط برنامه نویسی نرم افزار ۳۵20۷۵۲۷ ۳0۱/۵۲ "1۳" فضایی را برای انجام مکانیزه برخی از 
کارها فراهم می‌کند که ساختار استفاده شده در */]۳([** شبیه ساختار دستوری زبان برنامه نویسی +/) است. در این پژوهش. با 
بازیابی بار به‌منظور دست‌یابی به اهدافی از جمله کاهش زاویه فاز ایستا و انرژی تامین نشده پرداخته شده است. مقدار توان بارهای 
شبکه( متغیرهای کنترلی) توسط الگوریتم بهینه سازی تولید شده و به محیط نرم افزار 1 تا ][9ع]1 انتقال می‌يابند. در این محیط 
این متغیرها باید در درون یک بردار تعریف شده و در زمان استفاده فراخوانی شوند. 
مراحل حل مسأله به شرح زیر است: 

۰ مرحله دوم: ثبت اطلاعات تولید شده در یک فایل واسط 

۰ مرحله سوم: خواندن اطلاعات در محیط ۱9۳۱۸ در این مرحله اطلاعات ثبت شده در محیط واسط توسط نرم افزار 

۱9۳۱۸ خوانده شده( مقادیر توان بارها) و در یک بردار از پیش تعبیه شده ذخیره می‌شود . 

٩‏ مرحله چهارم: انجام محاسبات پخش بار و تعیین اطلاعات مورد نیاز 
در این مرحله با توجه به محل نصب منابع و همچنین. ظرفیت آن هاء شبکه آماده انجام محاسبات پخش بار است. پس از انجام محاسبات 
پخش بار, میزان زاویه فاز ایستا و میزان آنرژی تامین نشده تعیین می‌شود. 

۲/۸1/۸3 مرحله پنجم: انتقال نتایج پخش بار به محیط نرم افزار‎ _ ٩ 


8 یک طرح بهینه جهت هماهنگی بازیابی بار و کاهش زاویه فاز ایستا با استفاده از الگوریتم (11,۳6[" 


ورود داده های مربوط به واحدهای تولیدی» اطلاعات 
شبکه» بارهای الکتریکی 


درنظر گرفتن حدود متغیرهای تصمیم گیری 


فرمول‌بندی مساله بازیابی باربا هدف کاهش زاویه فاز ایستا و 
انرژی تامین نشده 


حل مسئله برنامه ریزی چندهدفه با استفاده از 
الگوریتم بهینه سازی آموزش و یادگیری 


بررسی وجود انحراف از 
قیود وانجام اعمال اصلاحی 


چاپ نتایج 


شکل(۲): روندنمای پیشنهادی روند حل مسئله بازیابی بار 
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در این:مرحله اطلاعات توسط یک فایل واسط به محیظ نرم اقزار فلضب1۲ 0۸ متتقل شده و اجازه ادامه محاسیات به این نرم افزار داقه 
می‌شود. در طی روند حل مسئله بهینه سازی» هر زمان که نیاز به محاسبه تابع هدف باشد مراحل بالا برای تعیین مقدار تابع هدف طی 
می‌شود. شکل (۲) مراحل بیان شده در بالا را به طور کلی نشان می‌دهد. 


۳- شبیه سازی مساله پیشنهادی 

در این بخشء مسائل ارائه شده به‌همراه روش‌های حل داده شده در بخش قبل بر روی شبکه استاندارد ۳۹ باس اعمال می‌گردد و نتایج 
شبیه‌سازی آن‌ها مورد تجزیه و تحلیل قرار می‌گیرند. دو تابع هدف مورد نظر در این پژوهش کاهش زاویه فاز ایستا و کاهش انرژی تامین 
نشده می‌باشند. جهت بررسی وجود پا عدم وجود تضاد بین توابع هدف ابتدا به‌صورت مجزا هر کدام از توابع هدف «کاهش زاویه فاز 
ایستا (8۳۸) و کاهش انرژی تامین نشده(۳(/9) بهینه‌سازی گردیده و سپس نشان داده می‌شود که وجود تناقض بین این دو تابع هدف 
باعث می‌گردد که مسئله بازیابی بار به‌صورت چندهدفه و هم‌زمان مورد بررسی قرار گیرد. نتایج بهینه‌سازی تک‌هدفه نیز به‌عنوان مقادیر 
حدی در بهینه‌سازی چندهدفه مورد استقاده قرار می‌گیرد. بهینه‌سازی با برقراری اتصال بین دو نرم افزار دیگسایلنت(جهت پخش بار) و 
متلب(جهت بهینه‌سازی) صورت گرفته است. 


۱-۳- شبکه ۳٩۹‏ باس ,۲۳:۲۲ 

شکل(۴) دیاگرام تک خطی شبکه استاندارد ۳۹ باس 18۳ را نشان داده است. میزان تولید و بار کل شبکه ۳٩‏ باس به‌ترتیب براپر با 
۷ 6140 و ۷۲۷۷ 6097 می‌باشد. از آنجایی که در هر زير سیستم می‌بایست حداقل یک واحد خودراه انداز وجود داشته باشد با توجه 
به مرجع [۲۰] سیستم تحت مطالعه دارای ده واحد خودراه انداز می‌باشد که معیاری برای انتخاب تعداد زیر سیستم‌ها می‌باشد. بر اساس 
مرجع [۲۰] با قطع خطوط ارتباطی بین باس‌های شماره ۰۲۶-۲۵ ۱۸-۱۷ ۱۴-۴ و ۱۴-۱۳ دو جزیره مجزا تشکیل می‌شود. که در 
شکل (۴) نشان داده شده است. 


انجام بهینه سازی و خواندن اطلاعات و 


محاسبه متغیر های 
مسئله 


ی 


ارسال اطلاعات ارسال اطلاعات به 
به محیط واسط محیط و اسط ۳ 
شخایه اع قافن سا 1 
و( افت اطا ات : 
۳ در بافك َ 
توسط دیگسایافت س کت ۰ 
- مورد نیاز شبکه 3 


شکل(۳): فر آیند کلی حل مسئله توسط لینک نرم افزار دیگسایلنت و متلب 


3 فناوری‌های نوین در مهندسی برق و سیستم انرژی سبز سال اول. شماره ۱ بهار ۱۴۰۱ ِ" 


هادی حسین‌پوره محمدرضا اسماعیلی : امین خدابخشیان 
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شکل(۴): شبکه استاندارد ۳۹ باس 1۳:۳۳ [۲۰] 


1 ۶۱۷۲۵۲ ابو 


۲-۳ کاهش زاوبه فاز ایستا(8۳۸) و انرژی تامین نشده(۳۱۹) 

در این بخش. به‌منظور کاهش زاویه فاز ایستا و انرژی تامین نشده. توابع هدف (معادله (۶) و (۷)) به دو صورت تک‌هدفه و چندهدفه 
مورد تحلیل و بررسی قرار می‌گیرد. در مرحله اول طراحی تک‌هدفه برای هر یک از توابع هدف بیان شده در بخش دوم مقاله (معادله 
در مرحله دوم طراحی چندهدفه اجرا می‌گردد تا جواب‌های بهینه نهایی مساله (متغیرهای تصمیم‌گیری بهینه) شامل توان اکتیو و 
راکتیو بارها به‌دست آید. 


۳-۳- بهینه‌سازی توابع هدف به‌صورت تک‌هدفه 

در این بخش توابع هدف بیان شده (معادله (۶) و (۷)) در طول فرآیند بازیابی بار بصورت جداگانه بهینه‌سازی می‌شوند. در این راستا؛ با 
اجرای الگوریتم بهینه‌سازی مقادیر بهینه ماکزیمم زاویه فاز ایستا بین باس‌های شماره (۱۴-۴). (۰۱۴-۲ (۱۷-۱۸) و (۲۶-۲۵) 
وهمچنین انرژی تامین نشده در ۱۰ مرحله بارگذاری محاسبه می‌گردند. همان‌طور که بیان شد. جهت بهینه‌سازی توابع هدف (معادله 
(۶) و (0۷) از الگوریتم آموزش و یادگیری(11130) استفاده می‌شود. به‌منظور بررسی عملکرد این الگوریتم با الگوریتم‌های ابتکاری 
دیگر از جمله ۳۹0 [0۸۵۰]۲۱ [۲۲]و ۳۸ [۲۳] نیز مقایسه شده است که نتایج در ادامه آورده شده‌اند. با انجام بهینه‌سازی به‌صورت 
تک‌هدفه و در نظر گرفتن توابع هدف مربوطه. منحنی‌های همگرایی به‌تر تیب در شکل (۵) و شکل (۶) نشان داده شده‌اند. مقادیر 
مربوط به توابع هدف محاسبه و نتایج این طراحی در جدول (۱) داده شضده‌اند. اين مقادیر به‌عنوان مقادیر حدی مورد نیاز در طراحی 
چندهدفه مورد استفاده قرار می‌گيرند. 
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شکل(۵): مقایسه منحنی همگرابی کاهش زاوبه فاز ابستا با استفاده از روش‌های مختلف 
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شکل(۶): مقایسه منحنی همگرابی کاهش انرژی تامین نشده با استفاده از الگوریتم‌های مختلف 


جدول (۱): نتایج طراحی تک‌هدفه با استفاده از الگوریتم ۲1,80 


1۳15)0۷۲۷۷( 


(۳۸)682 نوع طراحی تک‌هدفه 
طراحی بر اساس تابع هدف ماکزیمم زاویه فاز ایستا(معادله ۶) 


۳۶۶۰ ۸.۰۱ طراحی بر اساس تابع هدف انرژی تامین نشده(معادله ۷) 


با انجام بهینه‌سازی در توابع هدف به‌صورت جداگانه. مقادیر بدست آمده برای هریک از توابع هدف در جدول (۱) آورده شده است. 
همان طور که در این جدول مشاهده می‌شود به ازای بهینه‌سازی به‌صورت تک‌هدفه و با هدف کاهش ماکزیمم زاویه فاز ایستا مقادیر 
٩۳۸-8۵‏ و 83711۷1۷۷۳-(۳ و به ازای بهینه‌سازی به‌صورت تک‌هدفه و با هدف کاهش انرژی تامین نشده مقادیر 
۳-۷ و ٩۳۸8/1068‏ بدست آمده است. نتایج جدول نشان می‌دهند. بهینه‌سازی هرکدام از توابع هدف تابع هدف 
دیگر را تضعیف می‌کند. یا به عبارتی دیگر با بهینه‌سازی هر کدام از توابع هدف. مقادیر بدست آمده در تابع هدف دیگر بدتر می‌شود. از 


این‌رو نتیجه می‌شود. این دو تابع هدف با هم در تعارض بوده و جهت دست یابی بهترین جواب بهینه برای تحقق هر دو هدف. باید 
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هادی حسین‌پوره محمدرضا اسماعیلی امین خدابخشیان 


مسئله را چندهدفه حل نمود. برخلاف روش‌های بهینه‌سازی تک‌هدفه. روش‌های بهینه‌سازی چندهدفه منجر به یافتن تنها یک جواب 
بهینه نمی‌ش‌وند. بلکه مجموعه‌ای از جواب‌ها را مشخص می کنند که تعادل مناسبی بین اهداف مختلف ایجاد می‌کنند. در ادامه. 
بهینه‌سازی چندهدفه توابع هدف مذ‌کور توضیح داده شده است. 


4-۳- بهینه‌سازی مساله پیشنهادی به‌صورت چندهدفه 


در این قسمت. طراحی چندهدفه مبتنی بر الگوریتم پیشنهادی (الگوربتم بهینه‌سازی آموزش و یادگیری) بر روی سیستم مورد نظر 
انجام شده است تا بهینه‌ترین مقادیر مربوط به متغیرهای تصمیمگیری(توان‌های اکتیو و راکتیو بارها) با بهینه‌سازی هم‌زمان دوتابع 
هدف مورد نظر (معادله (۶) و معادله (۷)) حاصل شود. منحنی همگرایی پرتو استخراج شده از این طراحی در شکل (۷) ترسیم شده 
است. در شکل ۷ ستاره ها (قرمز و آبی رنگ) جواب‌های غالب در آخرین تکرار حل مسئله می‌باشند. سازگارترین جواب در بین تمام 
جواب‌ها با در نظر گرفتن درجه اهمیت هر یک از توابع هدف توسط تصمیم گیرنده (بهره‌بردار سیستم) انتخاب می گردد. 

به‌منظور انتخاب یک جواب مطلوب. دو جواب متفاوت با در نظر گرفتن وزن‌های متفاوت برای توابع هدف مورد نظر به‌دست می‌آیند 
که نتایج آن در جدول (۲) آورده شده است. پارامترهای ۷۷1 و ۷۷2 به‌ترتیب نشان دهنده وزن‌های دوتابع هدف حداقل زاویه فاز ایستا 
و حداقل انرژی تامین نشده می‌باشند. معیار انتخاب این ضرایب با توجه به شرایط سیستم و اولویت‌های بهره برداری حاصل می‌گردد. 
در این مقاله فرض بر آن است که برای طراح هر دو تابع هدف از اولویت یکسانی برخوردار می‌باشد. با توجه به جدول اختلاف مقدار 
انر‌ی تامین نشده با در نظر گرفتن ضرایب با مقدار مشابه ۰,۵ (۴۸۳۴ مگاوات ساعت) با ضرایب 1۷2-07 , 1۷1-0/3 (۴۶۶۱ مگاوات 
ساعت) تنها حدود 3/5 درصد است. از طرفی دیگر حداقل مقدار زاویه فاز ایستا با در نظر گرفتن ضرایب مشابه ۰,۵ (۶درجه) در 


علاوه برکاهش ۲۴درصدی زاویه فاز ایستاء مقدار انرژی تامین نشده را در نزدیکی مقدار اولیه خود قرار می‌دهد. با در نظر گرفتن ضرایب 
۱۷2-3 ۳ ۱۷1-7 مقدار انرژی تامین نشده ۳ مکاوات ات نسبت به حالت دوم افزایش می‌یابد و زاویه فاز ایستا نیز حد‌ود 


8 یک طرح بهینه جهت هماهنگی بازیابی بار و کاهش زاویه فاز ایستا با استفاده از الگوریتم (11,۳6[" 


۰درصد کاهش می‌یابد. لازم به ذکر است که در اين مقاله وزن‌های در نظر گرفته شده برای توابع هدف موجود در طراحی به‌طور 
یکسان برابر با 5 (۱۷2<05-<۷۷1) می‌باشد. به‌عبارتی دو تابع هدف از اهمیت یکسانی برخوردار هستند. اما تصمیم گیرنده (بهره بردار 


سیستم) بسته به شرایط و اهمیت هر کدام از توابع» ممکن است هر ضریب دلخواه دیگری درنظر بگیرد. 


جدول (۲): نتایج مربوط به ضرایب متفاوت در طراحی چندهدفه 


٩۳۸ )068(‏ )۳۱5 جواب غالب پرتو 
۹/۸ ۶2۱ (۱۷۰-0۵/7 ۰ 0/3-,۱۷) 
۸/2 2۷:۲ (۵/6-<۷۷۵ , ۱20/4 ۷۷) 
5 ۰۸۳ (۷۷۰-0/5 , ۷۲۱-0/5) 
5/5 ۱« ۰۷۰۵/4 0/6 ۱۷) 
۳/۶ ۰۷۷ (۲۷۰-0/3 , ۲۲۱-۵/7) 


با در نظر گرفتن ضرایب یکسان برای همه توابع هدف. مطلوب‌ترین جواب در شکل (۷) به‌صورت ستاره آبی‌رنگ انتخاب می‌گردد که 
تعانج ان ین مین ودره قح ات ااب ساکا ترتع مصواب سا مایی اگم یه گنای انسامرمس شوه که به )زاف ۵ 
تابع هدف موجود در مسئله به‌طور هم‌زمان دارای مقادیر قابل قبول‌تری نسبت به دیگر جواب‌های غالب (ستاره‌های قرمز رنگ) باشند. 
در این جدول, مقدار بهینه (توان اکتیو و راکتیو) بارهای قابل بازیابی در مرحله دهم از بازیابی بار و مقادیر مربوط به توابع هدف نشان 
داده شده است. بر اساس جدول (۳()۳). مقادیر به‌دست آمده برای توابع هدف در طراحی چندهدفه مغایر با مقادیر مشابه در طراحی 


تک‌هدفه می‌باشد. زیرا؛ هر دو تابع هدف موجود در طراحی چندهدفه به‌طور هم‌زمان بهینه می‌گردند. 
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شکل(۷): منحنی همگرابی پر تو مبتنی بر 10 


جدول (۳): نتایج طراحی چندهدفه با استفاده از روش پیشنهادی (11/80) 


مقادیر بهینه توان اکتیو و راکتیو بارها به آزای بهترین جواب 


رچ0۱۵/7۸ 1 ۳۷۲۲ | شماره بای | بار 0۱۵/۸ ] (۳۵/۲۷ | شماره باس | بار 
9۹/65 24/5 ۷۳ ۳ 2/4 3122 ۳ ۱ 
58 308/6 ۳۴ ۲۲ 184 500 ۴ ۲ 
4/04 224 ۲۵ ۱۳ 8417 233/68 ۷ ۳ 
16/96 139 ۲۶ ۴ 15/1 522 ۸ ۴ 
7531 281 ۲۷ ۱۵ 10/68 0/91 ۱۲ ۵ 
27/5 20۹/22 ۳۸ ۶ 132/96 120 ۱۵ ۶ 
2/65 23/5 ۲۹ ۱۷ 2224 329 ۱۶ ۷ 
0/18 035 ۳۱ ۱۸ 216 1/38 ۱۸ ۸ 
249/5 1104 ۳۹ 1۹ 102/9 28( ۲۰ ۹ 
15/8 274 ۳۱ ۱۰ 


٩۳۸-6 82 
۳۱9-4834 ۷ 


*- نتیجه گیری 
در این پژوهش. طراحی فرآیند بهینه بازیابی بار با هدف کاهش زاویه فاز ایستا و انرژی تامین نشده و ب‌صورت یک مدل چندهدفه 
فرمول‌بندی گردید. برای بهینه‌سازی این مسئله اساسی از یک الگوریتم ابتکاری با نام لا[ استفاده شده است. از ان جایی که 
مسئله به‌صورت چندهدفه معرفی شده؛ برای یافتن جواب بهینه از روش بهینه پر تو استفاده شده است. زیرا در این روش دو تابع هدف 
به طور هم‌زمان بهینه می‌شوند. روش پیشنهادی بر روی شبکه ۳۹ باس 1۳۳ پیاده سازی و نتایج آن ضمن مقایسه با روش‌های دیگر 
مورد بررسی قرار گرفت. نتایج شبیه سازی حاکی از موثر بودن طرح پیشنهادی در خصوص در نظر گرفتن همزمان دو تابع هدف 
کاهش زاویه فاز ایستا و کاهش انرژی تامین نشده می‌باشد. با توجه به نتایج شبیه‌سازی استخراج شده از اعمال روش پیشنهادی بر 
روی سیستم مورد مطالعه. نتایج زیر ارائه می‌گردد: 
* درنظر گرفتن هم‌زمان دو تابع زاویه فاز ایستا و انرژی تامین نشده باعث می‌شود تا در مقادیر توان بالاتر و با زاویه فاز کمتری 
بتوان عمل سنکرون ساز را صورت داد. با توجه به نتایج مقدار زاویه فاز ایستا در برخی شرایط حتی تا ۴۰ درصد کاهش نشان 
۰ استفاده از الگوریتم میتنی بر آموزش 9 یادگیری (0ظر11) 9 با یک ساختار ساده در قلب طراحی مورد نظر حجم محاسباتی 


۸ 
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یک طرح بهینه جهت هماهنگی بازیابی بار و کاهش زاویه فاز ایستا با استفاده از الگوریتم 100" 
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در فرآیند بهینه‌سازی را کاهش و سرعت همگرایی جواب را افزایش می‌دهد. با مقایسه انجام شده با الگوریتم‌های دیگر افزایش 
سرعت همگرایی به میزان تقریبی ۱۰ درصد قابل مشاهده می باشد. 

استفاده از روش بهینه پرتو در حل مساله پیشنهادی باعث می‌گردد تا تمامی توابع هدف موجود در مساله به‌طور هم‌زمان بهینه 
گردند و در نهایت با انتخاب سازگارترین جواب مقادیر مطلوب‌تری برای توابع هدف به‌دست آید. در بدترین شرایط با در نظر 
گرفتن ضرایب ۵۰ درصدی برای توابع هدف یک جواب مطلوب با میزان زاویه فاز ایستا ۶ درجه و انرژی تامین نشده ۴۸۳۴ 
مگاوات ساعت حاصل گردید. 


نتایج شبیه‌سازی نشان می‌دهد که الگوریتم ابتکاری پیشنهادی عملکرد بهتری نسبت به رقبای خود در بهینه‌سازی مساله 


پیشنهادی ارائه نمود. 


مراجع 
اصم‌مه-ا باه ی 0عووها طمتاه۲ماوع۲ مورک ۵۵ج م۵ 1۳۵882۲۵ رمع .7 2۳00 ره۷۸ ۷۰ رعصه ۷۷ ,۲۱ رتا20 ,2 
:1 ,2018 ,4955-4960 .0۵ ولل) رل 0۱6۲۵66 و۵6( ۵9 ۵۲۱۲۵ و99 رقطام 5۱5 
۰ ۵۵ 10 
٩۹1716‏ ۵۵ 0 1۳۵۲۵۹۵۵095 لا رهام ممتاهماوه؟ موه ۳۵0۳۵ علصز۳. با قمع ۸0101 ,1۷۲ 
7۰ :01 ,1994 ,22-28 00۰ و10.1 ,9 ,۷۵۱ 
۲ مماوماهتاو مصاونلمممتامعو ۵۶ )ماهتا روزا1۵۳2 ۷۰ 20 روطصانا با بللة۷۷ .۶ ر01665 ]0۱011605-1 .ل 
5 ۴۵۷۷۵ ۶۵۵۰۱6 ۰ ,۵۵)00010821 0مووه-مصتما ون هتاهعوو خر :صمت۲2م۲۵5 مونره ۵۵9۷۵ [مأ لته 
۰ 0 :101 ,2104 ,381-390 .00 ,116 ۷۵۱۰ ,]76560۲6 
عصوا عصعام ممتاه)ماوع؟ مم‌اووو امن ۵۴ ممت19021122 ۹۱۷ ومیهام۱۷]۵1۵ یه 20 ره)۷۵ ,۷ .1 با بها۱0 .خر .خر 
11 ,2007 ,۰1322-11329 00 و22 ,۷۵ ر۵اه فسوص ۵ ۳۵۱۹۵0 رل رقطووه 0۳۱/۴۳۲۴۲ 
۰ --(1۱ 
٩۱2206 ۰‏ اعقا مت ممتاهمتاون ‏ ممتاجتماوم: ‏ 20م1؟ ‏ رصشا بو 4ص و۷ ۷ ینام .۲۲ 
۰ ((۱ ۲ آن)) 0زا و۵6 0۷ ۵۵۲۰۵۱6۵ ۱۱۵۱۵۵۵۵ ۷۵ 2014 7 رطمتاه ۲۵608۴ 
۰ 2 :101 ,2014 ,1075-7079 
5 ب 1 رطملاهماوع۲ موی ۵۱۷۵۲ ۵۴ مصتصصتا عمج باتاتطاجمامتن:ر مط ۳1 م۳۲۱2 ۲۰ 24 ب«0نانا .ظ .1 
۰ ۷ :101 ,2012 ,3-9 .00 10.1 ,28 ,۷۵۱ رک۵/ ٩‏ 20۷۵ 070 
مجح ماه اصمتعصهه مه مصتسامنتاده؟ تتافنطد واتلتنا متتاهعاه ۵۴ )مجمرص 0صج ممممنا۵8860؟ رلها1 ۷ ۱۷۰ 
:1 ,2000 ,196-207 00۰ ,10.2 و88 ,۷۵۱ یلیل ۱۵ 0 ۳۳۵۵۵۵0۵ رولوتاحصه تناها اقصمام امه 
زرررمر کز/۱ ۲ ۱۲۱۳ 
۶ 1۵5)0120108 1 ممتامعصومی مصاً رم1فوتصصعصه) ۵۶ فلورافصه متصهصون رصحقطظ ۱۷۰ قمع ت4صفمص‌همطکا ,۷( 
8 01۱۱۴۱۸/۵۳ 60 ۳661۳۱6۵ 0۷ ۵۵۲۵۷۵۵ ۵۱۵0۵ 7 رقطه 5۹ 00۵9۷۵۲ 1۳1016019۳860160 
۰ .109/00۳/۳ 10.1 :401 ,2004 1624-1624 .00 ,(040۳137513 ۱0۰ .021 1۳۳۴۲) 2004 
٩5/۵116‏ ۵۵۴ 01 ۵۱6۵60۲96 1۳ ر 1 رطمتامنلع5 ماعصض معقطظ عصت‌صهاه ما امهمتوورض مه رامات۸ .1۷۲ 
6۰ :01 ,2000 ,470-478 .00 و10.1 ,9 ,۷۵۱ 
بل ۲ رطما۵؟۲۵9)0 ورد ۵۵9۷۵ 10۲ ممتامن۲۵0 همه معقطمر هصتل‌مهاو؟ رقطصتو نهر 220 ۳۱2227118 .(1 
-10.1049/10 :001 ,1998 82-88 .00 145 ما۷۵ و۵ وه موه ۳۳۵60-0۲۱ 
7۰ و 
۲۵90۲2600 موه ما ممتامنالع۳ مامصج معقطامر مصت‌صفاه »10 0مطامجهر ب«فط ۵ رتالنا ۷۰ 200 ۷۵ ,۲۲ 
۵۵ > 2و ا۵ع۵0ر ۵۴ هل ۱۵۲۱۱۵۵۵0۲۵۵ رقصحعجط امتاصمی ه وه متآم01 1020 10601001218 
...10009 [/16 10.10 :001 ,2013 ,664-674 .00 ,53 ۷۵۱۰ 9(/5/07016 
۲ 10۲ ۲۵0۱۲۵010۲ ماوصه معقطام عصتل‌صهاو ما طمهمممج سول بامالنباع۳ ۲۴ 4ص رتوطازصق؟ یه رتطامافیک .۸ 
,۰۳۲۰۲ ,301-307 .00 0.4 و12 ۷۵۱۰ ,۲عصاهو 20اه 00 ۱۳۵۱۹۵۵ 0۵۵0 رم۲2ماو۲6 اعاونوو 
10/100۰ :001 
مور ۵۳۷/۵۲ 1 رمتامناع1۲ مامصه-معهدام مصت‌صعاه مط رم میامتصمعع ۸ مصتصع ۷ ۲۰ مه متاممردملتطمطو .8 
1۱ ,2004 ,377-286 .00 و3 .10 و33 ,۷۵۱ ورگ کا۵0۵0 و 20 رظ۲۵۹]۵1۵10 
۳۹۴/۵ ۱۳۱۹۱۰۱/۵ 
۷۵۴۲ 1 امه معقطام عصال‌صهای هصتمن0ع۲ ۲۵۲ امن معهااه۷ عصتونا" راماانباه۳ ب مه 1صعصصلممطک ,1۷۲ 
:1 ,2017 و9-16 .00 و146 ,۷۵۱ ,۲۵9۵۵۲0 تک ۵2۲ ۵01۳6 "یطمتاه‌ماوع. طماونره 
۰ 1010101 
4ص مصت0[1 مه :۲۳۸۲9 رمجم‌هصصهن)-ومآمهصه .0 280 ر2م12-۳۴۵۲عصا مه ۲ راهبتباوو-ع) ۲۱۵۲ 5 بن) رقطمخ .۲ 
,۰ 079 ک ۷۷۱/۵۷ 07 ,۱607۷0۲15 ۷۵۲ ۵ج طا 0)قلناحطزو 
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]14[ 


۲ هر 1۱ دظ9 ۵۶ ممتامنالم؟ ۲۵۲ ۱۵۷1۵۵ 0۶۲۱۲۳۳۵ ۱79۵ معط عصتاهعتاوه1۳۳۷ ,اهنا .8 20 ۳۱۸22711۵ .(1 
:1 ,2018 ,479-491 .00 ,99 ,۷۵۱ و ۵6۳8 :(00010) و۵ ۵۴ مایا 0116 ول رصه‌اورد 
--10:1007/540031-018-0338 

۰ ۶0۱۱6۵071 ۱۷۵/۲۰۵۱۸۰۲ 1۵/۵ رولو همه موه ۵۷۵ م۰۱۷۵ بطاهتمها :1 فصه تتهطام .۲ ,1 
مط ‏ وعصنلع ۳۳۵۵۵‏ رهام ممتاهماوه؟ وروی ۵۵0۷۵ ۶۵۲ مه۲0 مره 1000812100 هر رثن ۰ 220 0۵10۲ .۲ 
47 ۵ 10 :01 ,2017 ,1234-1252 .00 ,20.7 ,105 ۷۵۱۰ 1۳۴۳۲۳ 

ما و٩‏ ۳۳۵9۷۵۲ رمصعاه: .ظ مضه رصقصطومم۳۱ یه نب رتلتقصو بک 3 مصهتطاوطل0202مط ,هر 
0۵۵6 0۳06012 ۱۵80072016 مصاتمل0زجصمن. فمصماوووماناه ۵ صمتاهمافع: افالمتهم تم تصممممه‌طوم 
تم ,2019, 1-19 .9.0 و4 .0ظ و30 ,۷۵۱ قاری و۱۵ 20066 ۵7 1۳۵5۵00 ۱۱6۲۱۱۵۵0۳9۵ 
۰( --(11 

۶ م0651 60امص0تزهوم ۵۷ مر رمصفتگ ۰ 20 رل‌صهصطونم۲۱۵ یه یک وطهتتافطل02021مطک یه متاتفحصوط .۴ ,1۷۲ 
مصماوری ممتفوتصوصهه ۵ واتاتمهتامد مومع قومممعور ممتاه؟0]وع1 1۵20 4ج محصصمطمو عصمتآمصمتاهعو 
:1 ,2018 ,23-37 ,00 ,102 ,۷۵۱ رکه و۵۲ وتو آ۵ع61۳ع ۵۴ ول ۱۱6۲۱۱۵۵0۳9۵ 
00092-17۰ [1. [/6 10.101 

0 م1617 ۷ ۳۰ 24 رطممصماهو ی رمطصت‌نامن ,۸ رفصتانه]۱۷ ۲۱۰ رقمعتزم ] اوطاصها‎ ۳۳۳2۲61018 ٩۱۷2۲۲۱ 
2001160 و16 .00 ماع ۴۵۲۵۲ 61:۵6 ۳۶ 2009 7 "یطمتاه0اوع۲ ماود ما‎ 2009, 11 
10(- 0۰ 

قداص عصتل‌مهاه 101 0مطامصه مصطتتمعله متامصمع ۸ رعصمطو ۷۰ 220 رصقناط ۷۰ رعصقطا2 .2 رانا .1 مطقتا ۷ .1 
11 ,2010 ,1-6 .00 ,۱۷۵۵17۵ 0۵7۵۳۵ فا لا "رجمتاج‌ماوه اوه بوسر طمتامتلع۳۵ موه 
۰ 10 

اتهای‌امفاها ورد مسمم لقتاون مامامم هد فعمتامو رفص ماماهمصم۳؟ روامله۳۲۱-20 ,۱۷ ,خر 
11 ,2014 و531-539 .00 و0.6ظ و9 ۷۵۱ وا ع ص1۳۵ 0۵۵۵0 1۳7 رط۲۵۹)۵1۵0 
۰ 20 -10.1049/161 
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صطتمعاه ممناهح‌نصنامن عصنتصععا فصه عصنطممع] 3 
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